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Abstract 

The charged lepton spectrum in semileptonic B — > X u lv decay near maximal lepton energy 
receives important corrections from subleading structure functions that are formally suppressed by 
powers of AQCD/ m b but are enhanced by numerical factors. We investigate the series of higher 
order terms which smear over a region of width AEp ~ Aqcd near the endpoint the contributions 
proportional to 5{Ei — nib/2) times (i) the matrix element of the chromomagnetic operator, and 
(ii) four-quark operators. These contribute to the total rate at the few percent level, but affect the 
endpoint region much more significantly. Implications for the determination of \V u b\ are discussed. 
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I. INTRODUCTION 



Extracting the CKM matrix element \V u b\ from the inclusive decay B — > X u lv is com- 
plicated due to the large background from the decay B — ► X c iv. A way to remove this 
background is to impose a lower cut on the charged lepton energy restricting it to be 
within a few hundred MeV of the kinematic endpoint. Unfortunately, in this region of phase 
space there are large corrections to the decay rate due to the b quark distribution function 
in the B meson, which introduces model dependence into the extraction of \V u b\. This un- 
certainty can be reduced by using the B — > X s 7 photon spectrum as an input |2], |3|, which 
was recently carried out by the CLEO collaboration M. 

The structure function which describes the lepton endpoint region in B — > X u iv and the 
photon endpoint region in B — > Xgj arises due to an infinite series of formally higher order 
terms in the operator product expansion (OPE) becoming leading in the endpoint region. 
Corrections to this description in the endpoint region are only suppressed by one power of 
(although nonperturbative effects in the total rate are suppressed by Aq CD /???|), 
but can be enhanced by other factors, as explained below. We concentrate on corrections 
which describe the smearing of certain terms in the OPE of order (Aqcd/^) 2 5(1 — y) 
and (Aqcd/^?)) 3 5(1 — y), where y = lEijra^. There are important terms at such orders, 
enhanced by numerical factors, which contribute to the B — > X u iv rate at the several percent 
level. Since these contributions are concentrated near the endpoint, their effects on the rate 
in a restricted region can be an order of magnitude larger. 

The B — > X u lv lepton spectrum, including dimension-5 operators |5| and neglecting 
perturbative corrections, is given by 



dr G 2 F ml\V uh \ 2 r 2 , Q Wl . 2A 2 
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where 



\ 1 = ^(B(v)\b v (iDfb v \B(v)), X 2 = ^(B(v)\b v ^a^G^b v \B(v)), (2) 

and b v denotes the b quark field in HQET. While Ai is not known precisely, the B* — B 
mass difference implies that A 2 = 0.12 GeV 2 . Of the terms in Eq. ([!]), the leading order (in 
A-QCD/ m b) structure function contains 2[0(1 — y) — Ai/(6m^) S'(l —y)+...]. The derivative 



of the same combination occurs in the B — > X s 7 photon spectrum || , given by 



oT = G 2 F m 5 b \V tb V*\ 2 aC 2 
dx 32 vr 4 



Ai — 9A 2 \ , , Ai + 3A 2 , , , Ai „ 



(3) 

where x = 2E 1 /m b . The effect of the leading order structure function on the semileptonic 
decay is reviewed in Sec. II. At subleading order, proportional to 5(1 — y) in Eq. (JJ) and 
to 5'(1 — x) in Eq. (|3]), possibly the most striking difference is between the terms involving 
A 2 , whose coefficient is 11/2 in Eq. ([I]) and 3/2 in Eq. (|3]). Because of the 11/2 factor, the 
A2 6(1 — y) term is important in the lepton endpoint region and it contributes about —5.5% 
to the total B — > X u £u rate. This contribution is investigated in Sec. III. 

At order Aq CD /t??|, there are dimension-6 four-quark operators in the OPE, 

V -A = (KYPlu) (uj^PlK), S -p = (b v P L u) (uP L b v ), (4) 

where Pl — (1 — 7s)/2. Compared to the lower dimensional operators discussed so far, these 

are enhanced by 167r 2 and also enter the B — > X u £u electron spectrum [0 as delta function 

contributions at the endpoint 

dr( 6 ) Gpml \V ub \ 2 2 . , 

= i27r f B ™<B(Bx- B 2 )8(l-y). (5) 

Here /b is the -B meson decay constant and B\ 2 parameterize the matrix elements of the 
operators in Eq. (Q) between B meson states, 

\(B\0 V - A \B) = ^ B U l -(B\Os- P \B) = ^ B 2 . (6) 

If one assumed factorization and used the vacuum saturation approximation then these 
contribution would vanish, since B\ = B 2 = 1 (0) for charged (neutral) B mesons. There is 
no reason to believe that factorization holds exactly, and deviations from it at the 10% level 
are quite possible [0]. With such an estimate (i.e., \B\ — B 2 \ = 0.1, and fp = 200 Me V), the 
four-quark operators contribute of order 3% to the total B — > X u tv rate. This contribution, 
which is absent from the B —>■ X s ^ rate, is investigated in Sec. IV. 

Section V contains our conclusions and the implications for the determination of \V u b\. 

II. THE LEADING ORDER STRUCTURE FUNCTION 

Near the endpoint region of the B — > X u iv lepton energy spectrum, the OPE is effectively 
an expansion in Aqcd/[^&(1 ~y)} instead of AqcdA 7 ^- Thus, within AEi ~ A QCD from the 



endpoint, there is an infinite series of terms of the form 

oT _ G 2 F \V ub \ 2 ml ~ A n 2ffM(1 A (?) 
dy~ 192vr3 to<ri 

which are equally important when smeared over Ay ~ Aqcd/^6 @- Here A n are given by 
the matrix elements 

\(B{v)\b v iD ai . ..iD an b v \B(v)) = A n v ai . ..v an + (other tensor strnctnres) , (8) 

and satisfy Aq — 1, A\ — 0, A2 = — Ai/3, etc. These terms can be formally resummed into 
a nonperturbative structure function for the b quark, 

f(k + ) = \{B{v)\b v 5(k + -iD-n) b v \B(v)) , (9) 

where n is a light-like vector satisfying n ■ v — 1 and n 2 = 0. The moments of f(k + ) satisfy 
/ dk + f{k + ) k™ = A n . This is convenient for model calculations because the effects of the 
distribution function can be included by replacing mj by ml = + k + , and integrating 
over k + , 



(10) 



where &T v /&Ei is the parton- level spectrum. In the region we consider in this paper, Eg. > 
2 GeV, taking the parton level spectrum as 2y 2 (3 — 2y) 8(1 — y) or just the leading twist part 
of it, 2#(1 — y), makes a negligibly small difference numerically, and therefore we will use 
the latter. (Perturbative order a s corrections are also neglected.) 

For purposes of illustration, we will use a simple model for the structure function which 
has three parameters, 

1 n ab k 

f(k + ) = x fp) (1 " x)ab ~ l e ~ a(1 ~ x) e{1 ~ x) ' X = T- 

The parameterization of Ref. || corresponds to b — > 1 and a — > a + 1. The b parameter 
is introduced because the first moment of the subleading structure function that occurs in 
Sec. IV need not vanish, contrary to the leading order structure function in Eq. (|9|) and the 
subleading structure function that follows from the one-gluon matrix element discussed in 
Sec. III. The first three moments of f(k + ) are given by 

J dk + f(k + ) = 1 , J dk + f(k + ) k+ = A{l-b), 

fdk + f(k + )k 2 + = A 2 \(l-b) 2 + -}. (12) 



Lepton energy- 


Fraction of events 


in the endpoint region 




intervals 


No smearing 


Ai = -0.1 GeV 2 


A 1 = -0.3GeV 2 Ai 


= -0.5 GeV 2 


2.0 GeV < Ei 


32.6% 


32.6% 


32.9% 


33.4% 


2.1 GeV < E e 


24.3% 


24.3% 


24.8% 


25.7% 


2.2 GeV < Ei 


15.9% 


16.1% 


17.2% 


18.3% 


2.3 GeV < E £ 


7.5% 


8.4% 


10.2% 


11.6% 


2.4 GeV < E t 


n/a 


2.6% 


4.6% 


5.9% 



TABLE I: Fractional rate in the endpoint region without and with resummation. For the structure 



function we used Eq. (11) with b = 1, A = 0.5 GeV, and Ai as shown in the Table. 



For the leading order structure function we will use Eq. ( |TT| ) with b = 1, A = 0.5 GeV, and 
\ 1 = -3A 2 /a = -0.3 GeV 2 , which are in the ballpark of the recent CLEO determinations [BJ . 

The effect of resumming the leading singularities is illustrated in Table |. For the struc- 
ture function we used Eq. fll!]) with b — 1, A = 0.5 GeV, and three different values of 
Ai(= — 3A 2 /a). This Table illustrates that the effect of the resummation becomes small as 
the lower limit of the lepton energy interval is lowered towards 2 GeV. Note that we have 
neglected perturbative corrections (this is probably the main reason why the numbers in the 
Ai = —0.3 GeV 2 column differ from those in Ref. Q). The value of the A parameter also 
affects the importance of resumming the singular terms. Of course the precise determination 
of A (i.e., the b quark mass) is crucial for any extraction of \V u b\ from inclusive B decays. 



III. SMEARING OF THE A 2 3(1 - y) CONTRIBUTIONS 

In this section we investigate terms proportional to A2#(l — y) in the B — > X u lv lepton 
spectrum (or to — x) in the B — > X s 7 photon spectrum), which are subleading in the 

twist expansion, and therefore do not necessarily drop out from the relation between the 
two spectra. Such corrections originate from three sources (for details see fTOfl): 

(i) corrections to the vanishing of the forward scattering matrix element of any dimension- 
4 operator, (B(v) \ b v iD a b v \B(v)) = 0, due to higher order terms in the HQET Lagrangian; 

(ii) corrections to the b(x) = e~ inibV ' x b v (x) relation between the b quark field in QCD and 
in HQET; 
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FIG. 1: The one-gluon matrix element. 



(iii) corrections to the leading contribution of the one-gluon matrix element in Fig. [I]. 

Correspondingly, there are three series of higher dimensional terms in the OPE that can 
be resummed into three subleading structure functions that describe the smearing of these 
terms in regions of widths AE? ~ Aqcd near the endpoint. These subleading structure 
functions that smear the (i), (ii), and (iii) terms discussed above are denoted in Ref. JTTJ] by 
t(u), hi(u), and /i2(^i, ^2), respectively. Since none of them are known, we will focus on 
possible effects of only one of them, which has the largest coefficient. 

Of these subleading structure functions, the one related to (i) is universal and may be 
absorbed into the redefinition of the leading order structure function. Therefore, this func- 
tion drops out from the relation between the B — > X u £u lepton spectrum and the B — > Xg'-f 
photon spectrum. To gain insight into the relative importance of the other two structure 
functions, note that the coefficient of the (Xz/ml) 5'(1 — x) term in Eq. (H) arises from the 
(i)-(iii) terms respectively as 3/2 = (3 — 2 + 2)/2 [[□]], whereas the corresponding decom- 
position of the coefficient of the (A 2 /mg) 5(1 - y) term in Eq. (0) is 11/2 = (3 + 2 + 6)/2. 
Therefore, we will focus on the subleading structure function that smears the one-gluon 
matrix element (iii). 

The one-gluon matrix element in Fig. [I] is given by the imaginary part of the B meson 
matrix element of 



bin 



1uP L b, (13) 



(rribV + k — q) 2 (rribV + k' — q) 2 

where e is the gluon polarization vector. The gluon field strength arises from the part anti- 
symmetric in the gluon momentum (p = k — k') and polarization, p^T A e Au * — > —{i/2)G^ u . 
The series of most important terms in the OPE in the endpoint region is obtained by ex- 
panding in k and k', and keeping the part where in each additional order in AqcD/^b the 
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Lepton energy 
intervals 


Fractional correction, 5T /V 
a = 2.5 a = l a = 10 


2.0 GeV < Ei 


-9% 


-9% 


-10% 


2.1 GeV < E t 


-12% 


-11% 


-13% 


2.2 GeV < Ei 


-16% 


-15% 


-18% 


2.3 GeV < Ei 


-24% 


-23% 


-27% 


2.4 GeV < E t 


-38% 


-42% 


-33% 



TABLE II: Fractional corrections to the rate in the endpoint region due to the A2 (5(1 — y) term 
from the one-gluon matrix element. The subleading structure function is modelled by Eq. ( |i~T| ) with 
b = 1, A = 0.5 GeV, and the different values of o shown in the Table. For the leading structure 
function a = 2.5 is held fixed (corresponding to Ai = —0.3 GeV 2 ), and 6 = 1 and A = 0.5 GeV. 



contribution to the rate becomes more singular |2], These terms enter proportional to 
matrix elements where each additional covariant derivative gives v a tensor structure. In 
analogy with A 2 in Eq. (0), we define 



]T -(B(v)\b v iD ai ...iD c 

m=0 ^ 



iD 



Otm+l 



iD an b v \B(v)) 



3(n + 1) A 2 X n v ai . . . v an + (other tensor structures) 



(14) 



The normalization is chosen such that Xq = 1, and Xi = because of the heavy quark equa- 
tion of motion. The higher order matrix elements are unknown. The resulting contribution 
to the B — > X u £i> lepton spectrum is 

dr Gl\V, lh \ 2 r> 3 



dy 192 7r 3 £T mi m 



00 y 

6A 2 E 



(15) 



As for the leading order structure function contribution, the effects of Eq. (|T3|) can be 
modelled by smearing the n = term with a function f(k + ), as in Eq. (|T0|), where the 
moments of f(k + ) must satisfy / dk + f(k+) = X n . 

In Table || we give numerical estimates of the effect of these corrections on the B — > X u lv 
rate in the endpoint region. We use for the leading order structure function the model in 
Eq. (1TT| ) with b — 1, A = 0.5 GeV and Ai = —0.3 GeV 2 . Since the leading structure function 
smears 6(1 — y), the sensitivity to the precise values of its parameters is relatively smaller 
than in the case of the subleading structure functions that smear contributions proportional 
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FIG. 2: Time ordered product that yields four-quark operators and higher dimension terms. 

to 5(1 — y). For the subleading structure function we use Eq. ( |TT|) with b — 1, A = 0.5 GeV 
and with different values of the a parameter. Using the B — > X s 7 photon spectrum as an 
input only cancels a third of these corrections, and furthermore, there is another subleading 
structure function that contributes comparably and with the same sign as the one shown 
in Table |TJ. The sign of these corrections is negative in all entries in Table [TI|, because the 
subleading structure function is assumed to be positive for all values of k + and because the 
contribution of this correction to the total rate is negative. 



IV. SMEARING OF THE FOUR-QUARK OPERATORS' CONTRIBUTION 

There are also important contribution to the lepton endpoint region from the four-quark 
operators in Eq. (|4]). Compared to the b quark decay rate, these contributions in the lepton 
endpoint region are suppressed by Aq CD /m^ but are enhanced by 16tt 2 (their suppression 
is Aq CD /77t| in the total rate). There is again an infinite series of higher dimensional terms 
in the OPE which are equally important in the region of lepton energies within order Aqcd 
from the endpoint. These are obtained from the amplitude in Fig. |2|, which is proportional 
to 

f 'P L i{m b i + z)i»Pd 



[u lv P L b) (16) 



{m b v + k — pe) 2 

The series of most singular terms which are equally important in the endpoint region of 
the lepton spectrum is obtained by expanding in k and taking the imaginary part, 

dF 00 2 n 

— oc £ - y) (<r - v^r) (v - p e r -..(«- ViY n 

dy ^ n\ ml 

x(B(v)\(b v ^P L u)iD ai . ..tD an (u lu P L b v )\B(v)} . (17) 
Here ftp = pe/mb, and the (g^ — v^v u ) term arises because in the kinematic region we are 



S 



Lepton energy 


Fractional correction 


\ST\/T 


intervals 


b = 1/2 


6 = 1 


b = 3/2 


2.0 GeV < Et 


9% 


9% 


8% 


2.1 GeV < E t 


12% 


12% 


10% 


2.2 GeV < E t 


17% 


16% 


12% 


2.3 GeV < E t 


28% 


23% 


15% 


2.4 GeV < E t 


57% 


37% 


18% 



TABLE III: Fractional corrections to the rate in the endpoint region due to the smeared four-quark 
operators. The first moment of this subleading structure function need not vanish, so it is modelled 



by Eq. (11) with A = 0.5 GeV, a = 2.5 and the different values of b shown in the Table. For the 
leading structure function we used b = 1, A = 0.5 GeV, a = 2.5 (corresponding to Ai = —0.3 GeV 2 ). 

considering v M dotted into the lepton tensor gives a subleading contribution. Since v — pe 
defines a light-cone direction, (v — pi) 2 yields a suppression by at least one power of Aqcd/^6- 
Therefore, we may write 



-(B(v)\(b v ^P L u) iD ai . ..iD an (u^P L b v )\B(v)) QT - v^v" 
= Y n v ai . . . v an + (other tensor structures) , 



(18) 



where the other tensor structures yield subleading contributions. Collecting all factors, the 
contribution to the lepton spectrum is 
dr _ G 2 F \V uh \ 2 ml f ^ 



dy 



3tt 

2G 2 F \V ub \ 2 m 2 b ~ Y n 



3tt 



(^ _ v ^ ( B (v)\(b v ^P L u) 5(l-y-iD- n) (u lu P L b v )\B(v)) 
Y^^(l-y). (19) 



n=0 " b b 

The first term in the expansion of the delta function reproduces Voloshin's result in Eq. (^j), 
since Y = {{B\O v _ A \B) - (B\0 S - P \B))/2 = {B l - B 2 )f 2 m B /8. 

An important difference between the structure function for the four-quark operators and 
those discussed explicitly in Sees. II and III is that in the present case the first moment 



of f(k+) need not vanish. The importance of this difference is illustrated in Table III. 
To obtain these estimates we held the leading order structure function fixed, as discussed 
previously. For the subleading structure function we use Eq. ( jTT|) with A = 0.5 GeV and 
a = 2.5 held fixed, while varying the b parameter which determines the first moment. Since 



such contributions do not occur in radiative B decay, using the B — > X s 7 photon spectrum 
as an input does not reduce these corrections. The sign of this correction is undetermined, 
and the magnitude is directly proportional to the assumed violation of factorization (we used 
\Bi — B 2 \ = 0.1 and Jb = 200 Me V). It may be possible to address the size of factorization 
violation using lattice QCD [Tj|. 

V. CONCLUSIONS 

We investigated subleading structure function contributions to B — > X u lv that are im- 
portant in the region Ef > 2 GeV, and could each alter significantly the relation between 
the lepton spectrum and the B — > X s ^ photon spectrum used for the extraction of \V u b\. 
We concentrated on two subleading structure functions. 

The first series of corrections, considered in Sec. Ill, is related to the smearing of the part 
of the ^2/77^) <5(1 — y) contributions to the B — > X u iv lepton spectrum originating from 
the one-gluon matrix element in Fig. 1. This is formally an order AqcdA^ correction in 
the endpoint region, however it enters with an enhanced coefficient, 6. Its contribution to 
the total B — > X u lv rate is known to be —3%, and estimates of how it affects the lepton 
spectrum integrated over different intervals near the endpoint are shown in Table f|. 

The second series of corrections, considered in Sec. IV, is related to the smearing of 
four-quark operator contributions. This is formally an order AQ CD /m^ correction in the 
endpoint region, however it is enhanced by 16n 2 . Contrary to the previous case, the sign 
of this contribution is unknown, and its magnitude is also uncertain. It could be a ~ 3% 
contribution to the total B — > X u lv rate, and estimates of how it affects the lepton spectrum 
integrated over different intervals near the endpoint are shown in Table |TTT| . 

These corrections could each affect the value of \V u b\ extracted from the lepton endpoint 
analysis above the 10% level, even when it is combined with the measurement of the B — > 
X s 7 photon spectrum. While varying the lower end of the energy interval between 2 GeV 
and 2.4 GeV, as done in the recent CLEO analysis, provides some indication that these 
corrections are probably not much larger than our estimates, the uncertainty is reliably 
reduced only if the lower limit of the lepton energy interval is near 2 GeV, in which case 
the lepton energy cut no longer eliminates the charm background. However, if the B — > 
X u iv lepton spectrum can be measured accurately down to 2 GeV, then it is likely that a 
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resummation of singular terms in the OPE in the endpoint region is not necessary. In that 
case, knowledge of the B — > Xg'j photon spectrum is not needed 1 and \ V u b\ can be extracted 
without a resummation of the leading singular terms. In the Eg > 2.2 GeV region, where 
the resummation of the leading singularities in the OPE is important, so is the uncertainty 
related to the subleading structure functions discussed in this paper. With more precise 
B — > X u lv and B — > Xg'-f data in the region above 2.2 GeV, it may be possible to observe 
these corrections. 
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